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ABSTRACT 



An elementary method of determination of the character of the hot phase transition 
in 4d four-fermion NJL-type models is applied to non-supersymmetric and supersym- 
metric versions of simple NJL model. We find that in the non-susy case the transition 
is usually of the second order. It is weakly first order only in the region of parameters 
which correspond to fermion masses comparable to the cut-off. In the supersymmetric 
case both kinds of phase transitions are possible. For sufficiently strong coupling and 
sufficiently large susy-breaking scale the transition is always of the first order. 
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1. Introduction. Recently a lot of attention has been devoted to the problem of phase 
transitions in a variety of particle physics models. A sound motivation for this activity 
is given by cosmology. In particular, it is believed that first order phase transitions 
imply several interesting and potentially observable effects. One well known and widely 
considered example is the possibility of baryon number generation during first-order 
electroweak phase transition [1]. Other examples include gravitational wave production 
during bubble collisions [2] and modifications of relic dark matter particle abundancies 
[3]. 

As the symmetry breaking in the standard model is achieved by means of a scalar 
field with a nonzero vacuum expectation value, it is the temperature phase structure 
of the bosonic models, for example 0^ theory, which attracts at present attention of 
the public and is subject to a careful investigation. At the same time similar prob- 
lems in simplest fermionic models, like the Nambu-Jona Lasinio (NJL) one, remain 
rather poorly known in four dimensions. It should be stressed that the fermionic mod- 
els, although generically nonrenormalizable in 4d, are very useful as phenomenological 
effective Lagrangians. They facilitate quantitative studies of phenomena for which 
fundamental theories are not known and provide a tractable description of dynamical 
symmetry breaking. On more practical side, a specific interest in this type of theories 
may be justified by composite Higgs models of BHL-type [4], [5] and by recent revival 
of studies on NJL- type QCD condensates [6] . 

In this note we shall discuss the temperature phase structure of the NJL and super- 
NJL (SNJL) models to leading order in large N expansion. An elementary but effective 
algorithm is used to determine the type of hot phase transitions in those models. We 
study parameter space of the nonsupersymmetric and supersymmetric four-fermion 
theories with particular regard to regions interesting for phenomenological applications. 

The structure of the paper is as follows. In section 2 we describe our approach 
to hot phase transitions in NJL models, in section 3 we apply the algorithm to the 
nonsupersymmetric model which allows for a fully analytical treatment, and in section 
4 we elaborate on the more complex supersymmetric NJL case. 
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2. NJL models of fermion condensation at finite temperature. We set out to inves- 
tigate the 4d four-fermion models to the leading order in large N expansion. Let us 
recapitulate briefly main steps of the standard procedure which converts the original 
Lagrangians into the form suitable for expansion. The goal is achieved through 
the introduction of auxiliary scalar fields [7], [8]. The original Lagrangian of the NJL 
model then readsS 



,2 



= id^X+a^X+ + idf.X.a^'X. - B*B + gBX_X+ + gB*X+X_ (1) 

where A_|_, A_ are chiral Weyl components of A and B is an auxiliary scalar field, non- 
propagating at tree level. In order to apply the expansion to this lagrangian one 
considers N generations of fermions A and rescales B \/NB, g g/\/N. Inte- 
gration over fermions gives an effective action for the auxiliary field B. The vacuum 
expectation value of that field is related to the fermion condesate through the relation 
< AA >= ^ valid to the leading order in Hence, the object of interest to us is 

the potential for B. To determine this potential it is sufficient to consider real valued, 
homogeneous B. Denoting m = gB one gets the effective action in the form 



SnjlH =N Jd^x - 2tr log(-92 + m^)^ 



(2) 



where S is written in the Euclidean space. We use finite temperature imaginary time 
formalism due to Matsubara i.e. we compactify the euclidean time on a circle with 
circumference (3 = 1/T. The fermionic functional determinant in is calculated with 
anti-periodic boundary conditions in the compact direction. 

The supersymmetric generalization of NJL model was constructed in [9]. The La- 
grangian with auxiliary superfields has the form 



:SNJL = y't/^^[($+<I'+ + $-$-)(! -A2^^M) + $2<I'2] 



+ 



j (i^6'[$i$2 - g^i^+^-] + c.c. 



(3) 



"'^We use notation of the book [13]. 



All superfields in (^) are cliiral, ^± contain fermions X^,X^, $1.2 are auxiliary super- 
fields and A is a susy-breaking scale. The first component of $1, which we denote by 
Ai, is proportional to < AA > and m = gAi is, as before, the mass of the fermion A. 
Functional integration over all matter fields leads to the following effective action for 
the parameter m 

SsNjdm] = N (^J d^x^ -2tr \og{-d^ + m^) + 2 tr \og{-d^ + + m^)^ (4) 

Performing direct evaluation of (0) and one obtains finite temperature effective 
potentials V^tjl and Vsnjl respectively 

+ A*log(l + ^)}-l/^log(l + e-^vW) (5) 



and 



+ AMog(l + ^)-K + A2)[A2 

- (m^ + A^)log(l + ^^)]-A^log(l + ^^^^)} 



(3 J (27r)3 

4 f d^k , -/3Vm^+A^+fc^^ 



where the divergencies have been regularized by a simple momentum cut-off at the scale 
A. One should note at this point that all the formulae we have written till now are well 
defined for all values of parameters. As discussed in the literature ([10], [11], [12]) the 
form of Vnjl and Vsnjl as functions of the temperature determines the character of 
the phase transition between localized vacuum configurations. Positions of the extrema 
of the potential are given by the gap equation which is just the condition for vanishing 
of its first derivatives with respect to m. 

Our purpose is to determine the type of the phase transition as a function of the 
parameters of the theory. First of all, let us make the following observation. By direct 



inspection of eqs. and one concludes that the derivatives of both effective poten- 
tials with respect to m consist of two parts: the first part is temperature independent, 
the second temperature dependent. It is easy to check that the temperature indepen- 
dent part is an increasing function of m while the temperature dependent one is a 
decreasing function of m, with monotonically decreasing (increasing) their first deriva- 
tives (i.e. the second derivatives of the respective terms in the potential), respectively. 
Moreover, the second derivatives of these functions (third derivatives of the respective 
terms in the potential) never cross each other. All together, this implies that each gap 
equation, which is the summ of the two terms, has at most two different from zero 
roots. The greater of them always corresponds to a minimum because the potentials 
are increasing functions of m for large enough m. Hence these effective potentials have 
at most one nontrivial minimum, at any temperature. 

There is another observation which simplifies the reasoning -the effective potentials 
of both models under consideration are analytic functions of around m = for 
nonzero T and A. This is due to the additional mass gap which appears for fermions 
compactified on a circle with antiperiodic boundary conditions. Hence one can expand 
the finite temperature effective potentials in the vicinity of m = in the Taylor series 
in the variable m^. 

We define the critical temperature in the standard way, demanding ^ q^"'' \m=o = 
= ^^^|rrt=o- Of course, below Tc the ^^|m=o is negative and above Tc - positive. 
At this critical temperature the term in the expansion which is proportional to can 
be negative or positive.l In the former case, at T = Tc the potential has to have a 
minimum for m 7^ 0. The phase transition which may take place at T 7^ is of the 
first order. This follows from the fact that if one slightly increases the temperature 
above Tc, one gets g^^|m=o > i.e. a "hill" arises between the local minimum at 
m = and the global minimum which exists for m > due to \ m=o,T=Tc < 0- In 

the second case the phase transition takes place ai T = and is of the second order. 

^We disregard the case when it vanishes as the one which requires a fine tuning of parameters. 
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Indeed, for T > the potential has a minimum at m = and shghtly below Tc one 
has g^^2) |m=o < which implies that a "valley" of an infinitesimal depth arises next 
to m = 0. There cannot be any other non trivial minimum for m 7^ as this would 
contradict the observation that the potential can have only one nontrivial minimum. 

In conclusion, we can determine the order of the phase transition by studying the 
behaviour of the effective potential in the vicinity of the point m = 0. 

In the next two sections we shall analyse phase transitions in NJL and SNJL models 
according to the prescription presented in this section. 



3. Nonsupersymmetric NJL model at finite temperature. Let us warm-up with the 
simpler, nonsupersymmetric model. Since we are operating in the vicinity of m = 
and the effective potential is an analytic function at that point we can Taylor expand 
it in m^, as explained in section 2. 

V^jL{m,T) = VMjLim = 0,T) + m''{-^ + \- ' 



12/52 87r2' 



l + 2c/ + 21og(/32A2 
327r2 



+ m^— '- + o{m''(3') (7) 



where = 1/T, c/ = 27 - 3/2 - 2 logvr ^ -2.64. 

According to the program outlined in section 2, we find the critical temperature 
demanding that the m'^ term in (^) vanishes at Tc. We get 

3A2 



T,^ = — (1-c) 

277^ 



where c - 

Next, the order of the phase transition is determined by the sign of the m*^ term 
at T = Tc. The analysis of the different physical situations is straightforward thus we 
simply summarize the results. For c > 1 - there is no phase transition as there is no 
(chiral) symmetry breaking at T = 0. For 1 — 2|!.gCj+o.5 < c < 1 - the phase transition 
is of the second order and it take place at Tc. For < c < 1 — 2|!.gC^+o.5 ^ q 22 _ the 
phase transition is of the first order and it may take place at T > Tc. One should note 
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that the squared temperature of the phase transition is greater then A^e^-''"'""'^ ~ A^/10 
in this case. 

Looking at the T = gap equation one can see that the last scenario corresponds 
to fermion masses comparable to the cut-off A. This may be relevant for QCD, [6], 
but it is definitely not the range of fermion masses suitable for top-mode standard 
model of BHL [4], where to have the top mass of the order of 10^ GeV one needs 
-^j2 ~ 1 — 0(10^^^). The top-mode scenario falls then safely into the region of the 
second-order phase transition (to the leading order in large N), which confirms earlier 
results [14]. 



4-. Susy-NJL model at finite temperature. The construction of the supersymmetric 
extension of the simplest version of the N JL model has been discussed carefully in [9] . It 
turns out that in exactly supersymmetric theory fermion condensation is impossible for 
any choice of parameters. However, if one introduces a soft susy-breaking scalar mass 
A, then the fermion condensate may appear and chiral symmetry breaking becomes 
possible. Following the same procedure as in the non-susy case we obtain to the leading 
order in large N, Taylor expanding in m^, the supersymmetric counterpart of 

VsNjL{m,T,/\) - Vsnjl{Q,T,/\) = 

4,1 + 2c/ A2 



327r2 167r2(A2 + A2) 
' ■log(|^)-T^/o(^^A^)) + o(mW (9) 



167r2 °'A2 + A2' 47r2 
where = dyy"'{y'^ + x)'^^"^ /{exl){^/y^~+~x — 1). 

Following the procedure applied already in section 3 we find the equation for the critical 
temperature demanding that the rri^ term in formula (P) vanishes: 

T^2+-2-f-2 l0g(l + ^VA^) + 1^2^2{P'A') = (10) 

The resulting equation cannot be solved algebraically due to the presence of the func- 
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tion I2. In order to give a flavour of tlie dependence of Tc on parameters of the theory 
we can find a rough approximate solution noting that values of the function I2 range 
from 7r^/6 to at 00. We arrive at the formula 



on the product (3 A. The condition for vanishing of this term can be easily solved. We 
get that if < A/2.16 the phase transition is of the second order. For > A/2.16 
the phase transition is of the first order. The characteristic scale which determins the 
critical temperature is now the susy-breaking scale. 

The analytic results one can get using the above approximations can be summarized 
as follows. If c > 1 - there cannot be any phase transition - the symmetry remains 
unbroken. If 0.07 < c < 1 - for small susy-breaking scale, < A < Ai, again 
no phase transition is possible; for A > Ai there is second order phase transition. 
Its temperature is given by eq.(^). The boundary value of A i.e. Ai is given by 
^^^(T = 0, Ai) = (eq.(0)). If < c < 0.07 - for < A < Ai no phase transition 
is possible (Ai is as above); for Ai < A < A2 there is a second order phase transition; 
for A > A2 there is a first order phase transition. The value of A2 is given by the 
condition that both terms on the r.h.s. of @ vanish. 

The exact results of the numerical study are presented in the figure 1. All values 
in this figure are scaled by the cut-off A. The curve "c/l" represents the dependence 
of dl = Ai/A on c ("til" is the solution to the eq.([lO|) at T = 0). The curve "rf2" 
represents d2 = A2/A. Thus, above the curve "rfl" no phase transition is possible, 
between "rfl" and "(i2" only the second order phase transition occurs, below "(i2" the 
first order phase transition occurs. The remaining curves are solutions to the equation 
(10) for T/A = 0.1, 0,2, 0.25 respectively. These are contours of the fixed critical 
temperature in the plane (d,c). 

The phase transition in the super-top- mode scenario, [5], described in the leading 
order of the large N expansion of the SNJL model, can be of the first order as long 




(11) 



We also see that for A << A the term proportional to in eq.(|^) depends only 
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as one expects that the critical temperature for the electroweak breaking is close to 
200 GeV and the supersymmetry breaking scale A lies in the TeV range. There 
Tc > A/2.16 and in addition for fermion mass 0(200) GeV we get A2 < A. 

5. Summary. In summary, a simple method has been used to determine the character of 
the hot phase transition in 4d four-fermion NJL model, both in its non-supersymmetric 
and supersymmetric version. We conclude that in the non-susy case the transition may 
be of the first order, however in the region of parameters which correspond to fermion 
masses comparable to the cut-off. In the supersymmetric case also both kinds of phase 
transitions are possible. For sufficiently strong coupling and sufficiently large susy- 
brcaking scale the transition is always of the first order. We applied these results to 
models under current investigation. We have found that the super-NJL scenario of the 
top quark condensation can exhibit the first order phase transition in the physically 
interesting range of parameters. 

Authors thank J. Wosiek and K. Zalcwski for useful discussion. 
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Figure caption 



Figure 1. Exact results of the numerical study of the expression (9). Label d denotes 
the ratio A/A, c represents -^jj- The curve "cil" represents dependence of dl = Ai/A 
on c being the solution to the eq. (p!oD at T = 0. The curve "(i2" represents d2 = A2/A. 
The remaining curves are solutions of the equation (10) for T/A = 0.1, 0,2, 0.25 
respectively. These are contours of the fixed critical temperature in the plane (d,c). 
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